The 3D incompressible Euler with a passive scalar : a road to blow-up? 
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, The 3D incompressible Euler equations with a passive scalar 6 are considered in a smooth domain 

SlCl 3 with no-normal flow boundary conditions u ■ h\gn = 0. It is shown that smooth solutions 
blow up in a finite time if a null (zero) point develops in the vector B = Vq x V#, provided B has 
£SJ , no null points initially : u: = curl u is the vorticity and q — u: ■ V9 is a potential vorticity. The 

presence of the passive scalar concentration 6 is an essential component of this criterion in detecting 
the formation of a singularity. 
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I. INTRODUCTION 



■ It is known that the 3D incompressible Euler equations have an array of very weak solutions [IHII , but whether 
' a singularity develops from smooth initial conditions in a finite time has been a controversial open problem for a 
generation P4l6j . Most numerical experiments are performed on periodic boundary conditions : the review in (l7j 
cites more than twenty of these. However, the aim of this paper is to study the blow-up problem in the context of 
the evolution of divergence-free solutions of the Euler equations u(x,t) together with a passive scalar 8(x,t) 

— = -Vp, — = 0, — = d t + u-V, divw = 0, (1) 
Dt Dt Dt w 

in a smooth finite domain SJcK 3 with no- normal flow boundary conditions u ■ fi\gci = 0. The inclusion of 0, which 
could represent any passive tracer concentration [l8|-[20| , allows the vector V9 to interact with the fluid vorticity field 



■ uj = curl u which evolves according to 



= u- Vu. (2) 

Dt K ' 

Formally, it is easily shown that the equivalent of potential vorticity q = u> ■ V0 is also a passive quantity : see [22[ 
for a more general discussion of potential vorticity in the geophysical fluid dynamics context. To show this we write 
what has become known as Ertel's Theorem as [21j 



3. Dq fDu \ ( D6\ 

which is no more than a re-arrangement of terms after an application of the product rule. Clearly 

^ = 0. (4) 

A result of Kurgansky [23Tj25j (see also j26l - |28j ) can now be invoked for any two passive scalars whose gradients define 
the vector 

B = VqxW0, (5) 
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in which case B turns out to satisfy 

DB 

=B •Vtt, divB = 0, divw^O. (6) 

Dt ' ' w 

The _B-ficld is a cross-product of two normals to the material surfaces on which 9 and q ride and is thus tangent to 
the curve formed from the intersection of the two surfaces [2(| [2?J ■ This result is formally true for the gradient of 
any two passive scalars riding on a divergence-free Euler flow and is not dependent upon the definition of q used in 
this context, although the passivity of q is, of course, dependent on this. The key point here is that B contains the 
gradient of uj (in projection) and two gradients of 9. We propose to exploit the fact that the evolution of B in © 
takes the same form as that of the Euler vorticity field in © or a magnetic field in MHD [29|, |3(| . 



II. STATEMENT OF THE RESULT 



The Beale-Kato-Majda (BKM) theorem is a key result for the 3D Euler equations @. A subsequent modification 
was proved by Ponce [3l| in terms of the rate of strain matrix (deformation tensor) defined by S — i (Vit + Vit T ) : 

Theorem 1 There exists a global solution of the 3D Euler equations u £ C([0, oo]; H s ) (~l C 1 ([0, oo]; H 3 ^ 1 ) for s > 3 
if for every T > 

/ ||5(r)||£- (B s ) dr<oo. (7) 
Jo 



Conversely, if there exists a time T for which 

i-T 



I 

■Jo 



||«5(r)||i«,( H 3) dr = oo, (8) 



then limt_>.T ||5(i)||z,°°(R3) = oo. 

In the original BKM-result @ ||td||L°°(R 3 ) replaced ||<S||l°°(r3). The proofs in [1, Hl[ are valid for flow in all f2 = M 3 
but the techniques used in those papers, such as Fourier transforms and the Biot-Savart integral, are not readily 



adaptable to no-normal flow boundary conditions u ■ n|ao = 0. Ferrari 32[ circumvented this difficulty by adapting 
some ideas from the theory of linear elliptic systems to achieve these results for the boundary conditions u ■ n|an = 0. 
For our purposes, ([8]) is the key blow-up result for our finite domain Q. 

While 9 itself is a constant of the motion its gradient can easily be shown to satisfy (likewise for Vg) 

\\V6(t)\\ L oe m < ||V0(O)|| L oo ( n)cacp / ||5|| L o. ( n)*-. (9) 

Jo 

Hence J Q ||iS||z,=(n) dr controls these gradients as it does solutions of the Euler equations. The main questions revolve 
around the occurrence of null points (zeros) in Firstly, initial data for \B\ must be free of null points : then a null 
can potentially develop either by maxima or minima developing in 9 or q or if Vg and V# momentarily align at some 
point. £]IIII contains an example of simple initial data and a domain Q with no null points for \B\. In the following, 
t* is designated as the earliest time a null point occurs in \B\. It is, of course, possible that B could blow up earlier 
than t* by some other mechanism. 

Theorem 2 On a smooth domain S] C I 3 with boundary conditions u ■ n\on = 0, with initial data for which 
\B(x, 0)| > and \\B (x ,0) ||_L°°(n) < oo, if there exists a smooth solution of the 3D Euler equations in the in- 
terval [0, t*), then at the earliest time t* at which \B(x, t*)\ = 

t* 

||<5||i«.(n) dr = oo . (10) 

o 

Conversely, if J„ T ||5||z,°o(n) dr < oo in any interval [0, T] then \B(x, t)\ cannot develop a null point for t £ [0, T]. 
Proof : On the interval [0, t*) first take the scalar product of §5§ with B, where \B\ = B 



D(B 2 ) 
Dt 



B Vm B . (11) 
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Thus, dividing by B 2 = B ■ B and then multiplying by ln_B, which could take either sign, we obtain 



1 D\\nB 

2 DT 

Now multiply by | ln_B| 2 ' m_1 ' for 1 < m < oo 

1 D\lnB\ 2m 



2 



(In 5) [B-S-B] . (12) 



2m Dt 



l nj B|2(m-i)( ln £) (b-S-B) , (13) 



and then integrate over the volume, invoke the Divergence Theorem and the boundary conditions on D, and finally 
use Holder's inequality to obtain 



d - [ \lnB\ 2m dV < [ llnBl 2 '"- 1 \S\ dV 
Jn Jn 



2m dt J n 



2771-1 



Using the standard notation ||^||l!>(Q) = (Jq |-X"| p dV) P , (fl"4]l reduces to 

d 



< [ J \lnB\ Zm dVj (J \S\ Zm dV) . (14) 
>i/p 



dt 1 

which integrates to 



lnB|| L2 7T7 (n) < \\S\\ L , m{n) (15) 



\\lnB(t)\\ L2mm < ||lnS(0)|| L2m(o) + f \\S(r)\\ L2m{n) dr . (16) 

Jo 

Since is bounded we take the limit m — > oo to obtain 

|| InB(t)|| L c. (n) < ||lnS(0)|Uo O(n) + / ||5(r)|| L - (n) dr . (17) 

Jo 

Provided B has no zero in its initial data the log-singularity at |£?| =0 causes the left hand side to blow up at t* 
thereby forcing J* ||iS||l=°(J2) dr — > oo as t — > t* . 

Finally, it is immediately clear from (|17l) that if J Q T \\S\\L°°(n) dr remains finite in an interval t G [0, T] then no null 
can develop in B. ■ 



Remark : The scalar product within q and the subsequent vector product of the two gradients within B produce 
a rich set of possibilities when zeros form in V# and Vg, while ||a;||£oo(Q) simultaneously blows up. For instance, in 
the case when W9 = 0, while ||w||.L°°(n) certainly blows up at t* , it is not clear whether Vg blows up or not because 
of the scalar product within q. However, if it does then this is consistent with the equivalent of ^ for Vq. In the 
case when a null forms through a maximum or minimum in q, any simultaneous blow-up in q would obviously have 
to happen elsewhere in the domain other than the null point. Likewise, inequality (|17p is consistent with a blow-up 
in B which, if it occurred, would again have to occur elsewhere than the null point. 



III. AN EXAMPLE OF INITIAL DATA WITH NO NULL POINTS 



An important question is whether there is initial data such that \B(x, 0)| > 0. We proceed to find a simple example 
of a set of initial data it on a finite domain C M 3 from initial data on us and 9 such that |JB| > and ||S||loo < oo 
for the elliptic system 

curlit = u;, divtt = 0, u • n|an = . (18) 

The usual methods, such as the Biot-Savart integral, are hard to apply on this domain but for the elliptic system in 
(|18p. Ferrari (32l | has shown that for given a vector u>, the velocity field u can, in principle, always be constructed. In 
the next paragraph this construction is performed in an explicit example in which £1 will be determined later. 
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Take the example us = (1, 1, 1) T : we firstly observe that there is a velocity field v = (z, x, y) T which satisfies 
divv = and cuilv = (1,1, 1) T but we cannot be sure that it satisfies v ■ h = for any given domain O. Therefore 
it needs to be modified to satisfy the boundary conditions. To do this we introduce some potential such that 

u = v + W(j). (19) 

Note that curlu, = (1, 1, 1) T . To guarantee that (jTHJ) holds, <f> must satisfy the Neumann boundary value problem 

= -(z, x, yf-h, (20) 

on 

which always has a solution on any smooth domain tt. Thus we have been able to construct a velocity field u 
corresponding to u> = (1, 1, 1) T , that satisfies the boundary conditions. For simplicity, now choose 9 = \(x 2 + y 2 + z 2 ) 
(say) and calculate q, Vq and V9 

q = x + y + z, V9 = {x,y,zf, Vg = (1, 1, 1) T , (21) 

and then B 

B = {z~y, x- z,y-x) T . (22) 

Note that \B\ = only on the straight line x = y = z = t for t 6 R. Hence \B\ > on any smooth, finite domain il 
that avoids this line : this is enough to achieve our goal. 



A0 = O, 



dn 



IV. CONCLUSION 



These results raise curious questions regarding the nature of 3D Eulcr flow with a passive scalar. Physically 9 could 
represent, for instance, the concentration of a dye or a quantity of fine dust added to an Eulcr flow. As a passive 
quantity it would be appear to be innocuous but its presence introduces the gradient V# which interacts with u> and 
thereby introduces the second passive quantity q = uj ■ V0 into the dynamics. The key result is the stretching relation 
for B in ©, where B is simply a vector tangent to the curve that intersects the two material surfaces for 9 and q. 

The first null point in \B\ then drives J ||<S||l°°(£2) dr — > oo through the logarithmic singularity. The presence of is 
therefore essential to this mechanism. This raises the question whether this singularity is of a fundamentally different 
type than those that are thought to develop in bare 3D Euler flow with no passive scalar present? 

The proof of Theorem [2] shows that it is essential that \B\ starts with no null points, which rules out the use of 
periodic boundary conditions because B = x V9 has zeros for every value of t. Hence a comparison with the main 
body of numerical experiments is not possible [l2l - fl6j , although it would suggest that a numerical examination of the 
occurrence and nature of null points might be fruitful with the boundary conditions used in this paper. 

A further variation on this problem is that of the 3D Euler equations with buoyancy, which can be written in the 
following dimcnsionless form 

— + 9k = -Vp, — = 0, divw = 0. (23) 

Dt Dt v ' 

9, while still passive, is physically a dimcnsionless temperature and appears because the density has beenn taken to 
be proportional to 9 in the Boussinesq approximation. This changes equation ((2} to 

— + V9xk = uj-\7u. (24) 

The extra term V# x k makes no contribution to equation Q and so q remains passive. However, the BKM criterion 
for this system on a finite smooth domain would need re- working because of the presence of the k9 buoyancy term 
in (|23l). 
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